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Abstract 

We propose a new method, emission-line eclipse mapping, to map the velocity 
fields of an accretion disk in position space. Quiescent dwarf novae usually ex- 
hibit double-peaked emission line profiles because of disk rotation. Since a part of a 
disk having a different line-of-sight velocity is successively obscured by a companion 
in eclipsing systems, they show time-varying line profiles. We calculated the time 
changes of the emission-line profiles, assuming Keplerian rotation fields (v v oc r -1 / 2 
with r being the distance from the disk center) and an emissivity distribution of 
j ocr~ 3 / 2 . We, then, applied the usual eclipse mapping technique to the light curves 
at each of 12-24 wavelengths across the line center to map the region with the same 
line-of-sight velocity. The reconstructed images typically exhibit a 'two-eye' pattern 
for high line-of-sight velocities, and we can recover the relation, v v oc d~ x l 2 , on the 
assumption of an axisymmetric disk, where d is the separation between the two 'eyes'. 
We will be able to probe the Keplerian rotation law, the most fundamental assump- 
tion adopted in many disk models, by high-speed spectroscopic observations with 8-m 
class telescopes. 

Key words: accretion, accretion disks - - binaries: close, eclipsing - - novae, 
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1. Introduction 

It is widely recognized that accretion disks play principal roles in various objects, includ- 
ing protostars, cataclysmic variables (CVs), X-ray binaries (XBs), and active galactic nuclei. 
Among them, CVs are the most ideal objects for accretion-disk research for the following rea- 
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sons. They are composed of a mass-losing, Roche-lobe-filling late-type star and an accreting 
white dwarf, and their orbital periods typically range between 1.5 and 10 hr. Such timescales 
are very convenient for ground-based observations. They are numerous; in other words, there 
are many CVs in our neighborhood. Also, CVs are bright in optical bands and exhibit a variety 
of time variations, some of which are closely related to various accretion-disk processes. Unlike 
XBs, moreover, the irradiation flux from the central parts of the disk and from the central ob- 
jects is not very large, which makes accretion physics much simpler than otherwise. Hence, CVs 
are useful laboratories to investigate accretion-disk physics. However, the size of the accretion 
disks in CVs is typically ~ R & , which is too small to resolve with ordinary telescopes; that 
is, we always collect all of the photons emitted from the entire disk surface (except for some 
special cases of young stellar objects). To test the disk theory with observations, therefore, we 
are obliged to compare the results only in terms of the integrated values, which cannot uniquely 
constrain the spatial structure of accretion disks. 

There exists, however, a novel method to indirectly 'resolve' the surface brightness dis- 
tribution over the disk plane. That is eclipse mapping (Home 1985, 1993), which makes use 
of the eclipse of high-inclination binary systems. Since a part of the disk is successively ob- 
scured by the secondary star during an eclipse, the eclipsing light curves contain information 
regarding the brightness distribution over the disk plane. We can reconstruct a two-dimensional 
(2D) brightness distribution from one-dimensional eclipse light curves based on the concept of 
the maximum entropy. [As for the case of quasars, one may use microlens events by stars in 
intervening galaxies (e.g., Yonehara et al. 1998), but this method is irrelevant here.] Similarly 
important is the technique of Doppler tomography (Marsh, Home 1988). Since the line profiles 
reflect velocity fields within the binary systems, including those of the secondary star and the 
gas stream to the disk, we can map the line-emissivity distribution in two-dimensional velocity 
space by analyzing the line-profile variations during the course of the orbital motions. This 
method, however, provides no definite information regarding the velocity fields in position (real) 
space. 

We, here, consider a third method, by combining the above two methods, to map the 
velocity fields in position space over the disk plane. The goal is to verify whether the disk 
rotation velocity (v v ) simply obeys the Kepler's law, vk = \jGM\jr (with M 1 and r being the 
mass of a compact object and the distance from the central object, respectively). Moreover, we 
may be able to detect deviations from Keplerian rotation due to spiral shock patterns (Sawada 
et al. 1986a, b) and/or non-negligible magnetic fields. Thus, we can say that mapping the disk 
velocity fields in position space is one of the most fundamental tests of the basic disk models. 
Nevertheless, nobody has yet completed such a method, to our best knowledge, partly because 
for that purpose we need to detect changes of emission-line profiles during eclipse (lasting only 
several tens of minutes), which requires an enormous amount of photons, too much to collect 
with 4-m class telescopes. Recently built 8-m class telescopes can do this job, however. In this 
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paper, we describe the basic methodology and a simple test to probe the disk velocity fields by 
using the eclipse of emission-line profiles during the quiescence of dwarf novae. 

Note that an apparently similar method has been proposed by Vrielmann et al. (1995); 
however, no detailed description is available and here we present an independent, slightly 
different method of our own. In section 2, we give basic considerations and then present the 
adopted model assumptions and procedures. Next, we explain our calculations of model light 
curves, assuming Keplerian rotation velocity fields, and present the resultant one-dimensional 
velocity profile obtained by the proposed method in section 3. The final section is devoted to 
discussion. 

2. Model and Procedure 

2.1. Basic Considerations 

Dwarf novae (DNe) are a subclass of CVs and are unique in showing repetitive outbursts. 
During the quiescent states between outbursts, DNe usually exhibit strong, double-peaked 
emission lines due to Doppler effects caused by disk rotation. Therefore, each line profile 
contains information about the disk velocity fields, although we cannot spatially resolve them. 
In eclipsing systems, however, a part of the disk, which has a different line-of-sight velocity, 
is successively obscured by a companion star, giving rise to time variations of the emission 
line profiles. We then apply the usual eclipse mapping technique to the radiation flux at 
each of multiple wavelengths across the line center to obtain the emissivity distribution of the 
region moving with the corresponding line-of-sight velocity on the 2D disk plane. With the 
radial emissivity distribution being specified, we are able to obtain the rotational velocity as a 
function of radius on the assumption of circular motion. 

According to Home and Marsh (1986) and Marsh and Home (1988), who calculated 
and displayed the contours of the constant line-of-sight velocity (t>i. .s.) on the disk plane and 
the corresponding wavelengths (5\/\ = V\. Q , s Jc with A being the line center wavelength), 
the shapes of the contours are distinct for high and low line-of-sight velocities: The large-ui. o s . 
contours exhibit 'two-eye' patterns and their separation increases with a decreasing line-of-sight 
velocity. The small-Wi. o s . contours, on the other hand, show more complex patterns similar to 
the shape of a dipole field. For a reference, we display the contours of the constant line-of-sight 
velocity in Keplerian rotation viewed at the binary phase 0.0 for the case of an inclination angle 
of i = 80° in figure 1. The critical velocity discriminating these two different patterns is given 
by ±v ip (R ont )sini, with R ont being the radius of the outer edge of the disk, which corresponds 
to the wavelengths of the two (blue and red) peaks in the line profile. Since the large-t>i. os . parts 
are much easier to deal with, we, here, aim to map the emissivity distribution of those parts. 
In other words, we focus on an analysis of the wing parts of the double-horned line profiles. 
Then, as mentioned above, the relation between the line-of-sight velocity and the separation 
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Fig. 1. Contours of the constant line-of-sight velocity on the disk plane for the case of an inclination of 
i = 80° and a mass ratio of q = 0.5. This picture corresponds to the case of N = 12, where N is the number 
of wavelength bins. 

between the two 'eyes' (referred to as d) provides information regarding the disk rotational 
velocity fields on the axisymmetric assumption. Since i, Mi, and other binary parameters are 
not usually precisely determined, we do not pay much attention to the absolute values of t>i. . s . 
but are concerned with the relative relationship (e.g., a power-law index, if represented as a 
power-law) between u Lo>Si and d. 

2.2. Disk Model 

We calculated the line-profile variations during an eclipse for a simple case and tested our 
method to see to what extent we can derive information regarding velocity fields. We considered 
a semi-detached binary system composed of a primary star with a mass of M\ and a secondary 
star with a mass of M 2 . We assumed: (1) a geometrically thin accretion disk is formed around 
the primary, (2) disk materials are orbiting on the circular orbits around the primary with the 
Keplerian velocity, and (3) the secondary star fills its Roche lobe. In calculations, we set the 
binary separation, a, to be a unit of the length scale. We further assumed a line emissivity of 



(see Home, Saar 1991; Home 1994) between the radius of the inner edge of the disk, i? in = 0.05, 
and that of the outer edge, R ont = 0.5 in units of binary separation, a. (If we perform Doppler 
tomography at the same time, it is not necessary to assume j(r), since we can derive it on the 
assumption of the axisymmetric disk.) 
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As for the binary system parameters, we adopted a mass ratio of q = M2/M1 = 0.5, a 
binary separation of a = R & , and an inclination angle of % = 80°. In physical units, the Keplerian 
rotation velocity is 

VK = „ 10 00 ( ( ' ) km (2) 

V r VM ; V10 10 cm/ V ' 

That is, we needed to cover velocity ranges of v ~ (1000-3000) sini km s -1 for the radius range 
of r ~ 10 9 -10 10 cm and M\ ~ M & . Corresponding Doppler shifts are S\/Xq = vsini/c ~ (0.003- 
0.01) sini (with c being the speed of light). 

2.3. Procedure of the Emission- Line Eclipse Mapping 

We divided the emission line profile into N = 12 or 24 sections with a constant wavelength 
interval, AA/Ao = 0.014/iV. We calculated the continuum flux with an equivalent width (EW) of 
100A, and each line flux was normalized with this continuum flux. In the following, F cont = 1.0 
and inline (A) refers to the normalized line flux with respect to F cont . 

For each section, an eclipse light curve was constructed from the binary phase <p = —0.15 
to 4> = +0.15 with 51 phase bins. We summed up the emission of the two sections (i.e. on the 
blue and red sides) , which have the same absolute line-of-sight velocity, to increase the number 
of photons and normalized it with its maximum value. In some model, we only used either side 
of the emission profile, since that would be better when the line profile is highly asymmetric. 
Then, we made an eclipse map using the PRIDA code (Baptista, Steiner 1991, 1993) for each of 
N calculated light curves with different wavelengths. As for the default map, we, here, simply 
used an axisymmetric map (discussed later). We fixed the inclination angle to be % = 80°. 

3. Reconstructed Images and Velocity Fields 

3.1. Basic Model 

We first calculated the expected light curves based on the line-profile changes for the 
simple cases of N = 12 and 24; i.e., we supposed that 12 (or 24) light curves at 12 (24) 
different wavelengths are available. Figure 2 show the images of the eclipsed disk (left) and 
the calculated emission line profiles for the cases of N = 12 (middle) and for those of N = 24 
(right), respectively, for five different binary phases; from the top, 0= —0.12,-0.06,0.0,0.06, 
and 0.12, respectively. Here, the binary phase of = corresponds to the mid-eclipse phase. 
Note that the binary rotates counterclockwise. 

There are several features which deserve attention. First, most of the line profiles exhibit 
double peaks, and the wavelengths of the two peaks do not vary significantly (see, e.g., those 
in the first, third, and fifth rows). In contrast, the line wings are mostly missing in the third 
row because of the eclipse of the innermost part of the disk, which has large rotation velocities. 
Second, the line profiles before and after the eclipse of the central region (i.e., in the second 
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Fig. 2. Eclipsed disk images (left), emission line profiles [-Fiine(A)] for the case of 12 wavelength bins across 
the line profile (i.e. N =12; middle), and those of N =24 (right), respectively, for five different binary 
phases; from the top to bottom, <j) = —0.12,-0.06,0.00,0.06 and 0.12, respectively. The disk is seen by a 
viewer at the bottom of the page. 
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and fourth rows) still show two peaks at the same wavelengths as those in the other panels. In 
addition, one more weak peak is found in the right panel with finer wavelength bins (N = 24). 
Such subtle spectral-line changes are important for successful emission-line mapping, thus a 
huge number of photons are required. Obviously, double-horned line profiles are smoother in 
iV = 24 than those in iV = 12. Finally, line profiles at one phase (say, <fi = <fro) arid those at 
another phase with opposite sign (<fr = — 0o) are exactly symmetric with respect to the line 
center. Thus, the light curves at one wavelength (say, at Ao + AA) are exactly symmetric in 
time with respect to the minimum with those at another wavelength at Ao — AA. This may not 
be the case in realistic situations because of the finite thickness effects and possible warps in 
the disk. 

We next display in figure 3 the original maps (left), synthetic light curves (middle left), 
and the reconstructed maps obtained by using both of the red and blue sides of the line profiles 
(middle right), and those reconstructed by using only the red side (right), respectively. We used 
N = 12 wavelength bins for these mappings. Since the light variations of the red and blue sides 
are symmetric, as mentioned above, we simply summed up the two contributions in the light 
curves to acquire good photon statistics. The middle panels show a diversity of light curves 
at different wavelengths. The upper two panels representing the cases of high velocities show 
similar variations to those of the continuum, but with a shorter duration, since they have no 
complicated shape of the emission region. As we look downward, the corresponding velocities 
decrease and the images start to show a complex pattern. Accordingly, the light curves are 
changed in a complex fashion. Especially, the light curves at the minimum is nearly flat in the 
upper three rows (although the duration of the minimum decreases downward), which contrasts 
rather complex variations in other two rows. This feature suggests that the high-velocity parts 
of the line-profile variations seem to be easier to derive any concrete information. 

Eclipse maps are constructed with 21 x 21 grids. The constructed maps displayed in 
the middle right panels of the upper two rows show the characteristic 'two-eye' patterns, which 
recover the qualitative features of the original maps. In the lower three rows representing the 
cases of low line-of-sight velocities, in contrast, the constructed images (depicted in the middle 
right panels) hardly recover the original map. This is because we used an axisymmetric default 
image in the reconstruction, whereas the line emission profile of the original map is far from 
being axisymmetric. 

To fix such a problem, we only use the red side of the line profile and show the results in 
the right panels of the lower three rows in figure 3. Interestingly, the right panels can recover a 
part of the eclipsed region on one side more clearly. In the third row, especially, the right panel 
shows a clear eye-shaped pattern, which is not seen in the middle right one. This, therefore, 
suggests that we are able to obtain the information of velocity for the case of v i. . S; = 700-1050 
km s -1 even in N = 12. 
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Fig. 3. Original cmissivity maps (left panels), synthetic light curves (middle left), and the reconstructed 
maps (middle right) for the cases of N = 12. The right panels in the lower three rows are reconstructed 
maps derived solely from the red side of the line profiles. Each row of the panels represents, from the top 
to bottom, the line-of-sight velocities of V\. . s . = 1400-1750, 1050-1400, 700-1050, 350-700, and 0-350 km 
s _1 , respectively. In the light curves, the solid, dotted and dashed lines show those of the sum of the blue 
and red sides, blue, and red side, respectively. The disk is seen by a viewer at the bottom of the page. 
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Fig. 4. One-dimensional distributions of the original (dotted) and the reconstructed (solid) line emissivity 
profiles along the x-axis for ui.o.s. = 1400-1750 km s _1 (left panels) and 1050-1400 km s _1 (right panels), 
respectively. The flux of the eclipse maps is multiplied by 25. 

3.2. Line Emissivity Distributions 

To see line-emissivity distributions more quantitatively, we illustrate in figure 4 cross- 
sectional views of the emissivity profile along the x-axis. The emissivity profiles of the original 
map are illustrated by the dotted lines, while those of the reconstructed map are by the solid 
lines. Although the reconstructed one shows a broader distribution, the peak wavelengths do 
not differ much in both panels. This figure confirms that the spatial positions of the two 'eyes' 
in the reconstructed images are to good accuracy equal to those of the original ones. That is, 
the reconstruction is quite successful at least as long as data of ~ 50 observational runs over 
one eclipse are available with good photon statistics. 

3. 3. Models with Large Number of Grids 

In figure 5, we present only the eclipse maps which show the characteristic 'two-eye' 
pattern for the case of N = 24. The left panels are for 21 x 21 grids, as in figure 3, while the 
right panels show finer maps with 41 x 41 grids. The data for ui. . B . = 1575-1750 km s _1 are not 
shown, since no convergence is reached in the latter case. Comparing the left panels with the 
right ones, we notice that finer maps represent more circular images. However, the positions of 
two 'eyes' do not differ significantly. 

3.4- The Rotation Law 

Finally, we calculated, for each line-of-sight velocity, the distance from the primary star 
in the real space based on the separation between the two 'eyes' on the reconstructed image; 
the result is plotted in figure 6. As is clear, the reconstructed line-of-sight velocities are roughly 
on the straight, dotted line with a slope in proportion to d~ 1 / 2 . That is, our proposed method 
is feasible to reconstruct Keplerian velocity fields, as long as good photon statistics is achieved. 

To see how observational errors affect our results, we repeated the same procedure, but 
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Fig. 5. Reconstructed maps for the case of N = 24, left panels show models with 21 x 21 grids, while 
the right ones show those with 41 x 41 grids. The line-of-sight velocities are, from top to bottom, v\. . s . = 
1750-1925, 1400-1575, 1225-1400, and 875-1050 km s _1 , respectively. The disk is seen by a viewer at the 
bottom of the page. 
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Fig. 6. Line-of-sight velocities obtained by the image reconstruction technique as functions of the distance 
from the primary star. The results for the case of N = 12 are indicated by '+', while those of TV = 24 are 
by ' x ', respectively, both obtained with 21 x 21 grids. '*' also shows the results for the case of N = 24 and 
41 x 41 grids. For a comparison, we depict the Keplerian velocity relation, v oc r -1 / 2 , by the dotted line. 



added random errors in calculating the light curves at each wavelength. Namely, we calculate 
the observable 
following steps: 



the observable line flux, F^, from the theoretically modeled flux, F^ del according to the 



1. F™ d *\\) = F^(\) + F cont , 

2. F^(X) = F™ dcl (A) + 6(S/N)" 1 V^T^A), 
3- F^(X) = F^ t (X) - F cont . 

Here, e represents the Gaussian noise (with the mean being zero and the standard deviation 
being unity); the amount of random variations is given by the signal-to- noise ratio (S/N). Note 
that we set -F con t = 1-0, and F™ c dcl was calculated so as to give an equivalent width of 100A(see 
subsection 2.3). 

We calculated the cases of S/N = 10, 20, and 50 as representative examples. We found 
that in some cases convergence is not always guaranteed in the eclipse mapping process, because 
the light curve is not symmetric when random errors are added. Nevertheless, we expected good 
results, even if we solely consider the cases with successful reconstruction. Figure 7 presents 
the results of the relatively noisy cases of S/N = 10. The number of wavelength bins is N = 24. 
It is impressive to see that the Keplerian relation is roughly reproduced, although we used only 
those results which have the clear 'eye' pattern. These results demostrate that our proposed 
method is powerful even for the case of S/N = 10, as long as the reconstructed eclipse map 
displays the 'eye' pattern. 
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Fig. 7. Effect of noise in the case of ^(wavelength bins) = 24. In the case of S/N = 10, the left panel 
shows the result with 21 x 21 grids for the reconstruction and the right panel does the one with 41 x 41 
grids. For a comparison, we also depict the Keplerian velocity relation, v oc r -1 ' 2 , with the dotted line. 

4. Discussion 



We propose emission-line eclipse mapping to map the velocity fields in an accretion disk. 
We studied the feasibility of emission-line eclipse mapping using model simulations. Our re- 
sults show that those regions having high line-of-sight velocities can be precisely reconstructed, 
while the low line-of-sight regions are not, as long as an axisymmetric default image is used. 
These reconstructed maps of high line-of-sight velocities show the characteristic 'two-eye' pat- 
tern. The corresponding line-of-sight velocity recovers the originally adopted relation, which 
is proportional to d~ ' 5 . It means that this method is a potentially useful tool to probe the 
Keplerian disk rotation from observations of the line-profile variations. 

Although the Keplerian rotation is a fundamental assumption in constructing the so- 
called standard disk model, it is not always trivial from a theoretical point of view. For example, 
a rather hot accretion flow, called ADAF (advection-dominated accretion flow), predicts sub- 
Keplerian rotation (Kato et al. 1998); i.e., v v < vk- In addition, the velocity fields are affected by 
the presence of large-scale magnetic fields and/or spiral patterns, which will produce deviations 
from the exact Keplerian velocity of the order of 5v/vk~ va/vk ~ c s / Vk ~ H/ R ~ 1/30 (where va 
and c s are the Alfven velocity and the sound velocity, respectively, and H/ R is the aspect ratio 
of the disk). To detect such small effects we need to collect an enormous number of photons, 
which are unavailable at the present. Thus, we should await innovations of the observation 
technique. 

For successful reconstructions, we need rather frequent observations over the eclipse pe- 
riods. Only 8-m class telescopes, such as Subaru, can meet this requirement, which is underway. 
The next issues are to improve the current method so as to minimize the total observing times 
by modifying the default images and to perform observations studies (e.g., Vrielmann et al. 
1995). 
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